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Abstract 

A two-dimensional arrangement of toothpicks is constructed by the fol- 
lowing iterative procedure. At stage 1, place a single toothpick of length 1 
on a square grid, aligned with the y-axis. At each subsequent stage, for ev- 
ery exposed toothpick end, place a perpendicular toothpick centered at that 
end. The resulting structure has a fractal-like appearance. We will analyze 
the toothpick sequence, which gives the total number of toothpicks after n 
steps. We also study several related sequences that arise from enumerat- 
ing active cells in cellular automata. Some unusual recurrences appear: a 
typical example is that instead of the Fibonacci recurrence, which we may 
write as a(2 + i) = a{i) + a{i + 1), we aei n = 2'^ -\- i (0 < i < 2'^), and then 
a{n) = a{2'' +i) = 2a{i) + a{i + l). The corresponding generating functions 
look like nfe>o(-'^ + ~^ + 2x'^ ) and variations thereof. 

Keywords: cellular automata (CA), enumeration, Holladay-Ulam CA, 
Schrandt-Ulam CA, Ulam-Warburton CA, Rule 942, Sierpihski triangle 
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1. Introduction 

We start with an infinite sheet of graph paper and an infinite supply of line 
segments of length 1, called "toothpicks." At stage 1, we place a toothpick 
on the y-axis and centered at the origin. Each toothpick we place has two 
ends, and an end is said to be "exposed" if this point on the plane is neither 
the end nor the midpoint of any other toothpick. 

At each subsequent stage, for every exposed toothpick end, we place a 
toothpick centered at that end and perpendicular to that toothpick. The 
toothpicks placed at odd-numbered stages are therefore all parallel to the 
y-eods, while those placed at even-numbered stages are parallel to the x- 
axis. 

Fig. [l] shows the first ten stages of the evolution of the toothpick struc- 
ture and Figs. [2j[3]show the structure after respectively 53 and 64 stages. 
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Figure 1: First ten stages of the evolution of the toothpick structure. 
The numbers of toothpicks in the successive stages, r(l), . . . , T(10), are 
1,3,7,11,15,23,35,43,47, 55. 

Let t{n) {n > 1) denote the number of toothpicks added at the nth 
stage, with i(0) — 0, and let T{n) := Y^^=o ^(*) ^^'^ total number of 
toothpicks after n stages. The initial values of t{n) and T{n) are shown 
in Table [T] These two sequences respectively form entries A139251 and 
A139250 in IB]. 

The first question is, what are the numbers t{n) and T{n)7 We start by 
finding recurrences that they satisfy. As the number of stages grows, the 
array of toothpicks has a recursive, fractal-like structure, as suggested by 
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Figure 2: The toothpick structure after 53 stages (there are T(53) 
toothpicks) . 
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Figs. [T] [2] |3] (For a dramatic illustration of the fractal structure, see the 
movie linked to entry A139251 in |S|.) In order to analyze this structure, 
we consider a variant, the "corner" sequence, which starts from a half- 
toothpick protruding from one quadrant of the plane. In ^2]we establish 
a recurrence for the corner sequence (Theorem and in q3 we use this to 
find recurrences for t{n) and T{n) (Theorem [2p CoroUarjTjsl . Section ^ 
gives a similar recurrence for the number of squares and rectangles that are 
created in the toothpick structure at the nth stage, and S|5] gives a more 
precise description of the fractal-like behavior and discusses the asymptotic 
growth of T(n). 

The recurrences make it easy to compute a large number of values of 
t{n) and T{n), so in a sense the initial problem has now been solved. 

However, the toothpick structure is reminiscent of another, simpler, 
two-dimensional structure, the arrangement of square cells produced by the 
Ulam-Warburton cellular automaton (Ulam [17], Singmaster [13], Stanley 
and Chapman fTF, Wolfram [22l p. 928]). For this structure there is an 
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Figure 3: The toothpick structure after 64 stages (there are T(64) = 2731 
toothpicks) . 

exphcit formula for the number of ON cells at the nth stage, and a simple 
generating function for these numbers, as we will see in S}6) Theorem [e] 

This hint led us to look for a similar generating function and an explicit 
formula for the toothpick sequence. Our first attempt was a failure, but 
provided an surprising connection with the Sierpinski triangle, described in 

The generating functions for the toothpick sequence and for a number 
of related sequences have an interesting form: they can be written as 

x(a + /3x)[](l+72;2'-i+5x2'), (1) 

k>e 

for appropriate integers a,/3,7, ^, e. In ^we describe the relationship be- 
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Table 1: The toothpick sequences and T(n) for < ti < 49. 



tween such generating functions and recurrences for the underlying sequence 
(Theorem [7| . The generating function for the toothpick sequence is then 
estabhshed in Theorem [Sl 

Generating functions of the form nfe>i(l + 9kx'') have been used in 
combinatorics and number theory for a long time (for a survey see [5]), 
but generating functions of the form ([T]) may be new — at least, until the 
commencement of this work, there were essentially no examples among the 
170,000 entries in [5]. 

The following section, ^ gives a general method for obtaining explicit 
formulas from the generating functions (Theorem |9| , and the particular 
formula for the toothpick sequence is given in Theorem |10[ 

Both the toothpick structure and the Ulam-Warburton structure are 
examples of cellular automata defined on graphs, and we discuss this general 



framework in { 10 We have not been able to find much earlier work on the 
enumeration of active cells in cellular automata — the Stanley and Chapman 
American Mathematical Monthly problem |15j and the Singmaster article 
[IBj being exceptions. We would appreciate hearing of any references we 
have overlooked. 

Of course, two well-known examples show that one cannot hope to enu- 
merate the active states in arbitrary cellular automata: the one-dimensional 
cellular automaton defined by Wolfram's "Rule 30" [19], [20], [22] behaves 
chaotically, and the two-dimensional cellular automaton corresponding to 
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Conway's "Game of Life" [2] is a universal Turing machine. 

However, we were able to apply our techniques (with varying degrees of 
success) to a number of other cellular automata defined on graphs, and the 
last four sections discuss some of these. Section [Tl] discusses a structure 
built using T-shaped toothpicks. Sections T2|T4 discusses variations on the 
Ulam-Warburton cellular automaton. SectionTl2l deals with the "Maltese 
cross" or HoUaday-Ulam structure studied in [17], as well as some other 
structures mentioned in that paper. Section [T3| considers what happens if 
we change the rule for the Ulam-Warburton cellular automaton of S|6] so 
that a cell is turned ON if and only if one or four of its neighbors is ON. 
The final section {[ 14 ) discusses what happens if we change the definition 
of the Ulam-Warburton cellular automaton to allow all eight neighbors of 
a square to affect the next stage. There is another variation that could 
have been included here, in which the rule is that a cell changes state if 
exactly one of its four neighbors is ON. Again we have a formula for the 
number of ON cells after n generations — see entries A079315, A079317 in 
[8j. Many further examples of sequences based on generalized toothpick 
structures and cellular automata are listed in [14] . 



Notation. Our cellular automata are synchronous, and we normally use 
the symbol n to index the successive stages. Cells are either ON or OFF. 
In all the examples we consider here, once a cell is ON it stays ON. Lower 
case letters (e.g. a{n)) will denote the number of toothpicks added, or 
cells whose state is changed from OFF to ON, at the nth stage, and the 
corresponding upper case letters (e.g. A{n)) will denote the total number 
of toothpicks or ON cells after n stages (the partial sums of the a{n)). By 
the generating function for a sequence a{n) (say), we will always mean 
the ordinary generating function A{x) := ^'^^oO.{n)x'^. If A{x) is the 
generating function for a{n), then A{x)/{1 — x) is the generating function 
for A{n). 



Remarks. 1. A common dilemma in combinatorics is whether to index 
the first counting step with n = or n = 1. In this paper we have consis- 
tently started the enumerations with zero objects (toothpicks or ON cells) 
at stage 0, adding the initial object at stage 1. This seems natural, and is 
the indexing used for most of these sequences in [8]. On the other hand, 
this is responsible for the leading factor of x in the generating functions ([T]) , 
(151, (17 1, etc., and for the fact that in the recurrences ([2|, (|4]), etc., the 



exceptional cases occur at the beginning of each block of 2^^ terms, rather 
than at the end. If they had occurred at the ends of the blocks, the be- 
ginnings of all the blocks would have agreed, which would have made the 
triangular arrays such as that in Table [3] look rather nicer (compare Table 
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[T]). Probably there is no perfect solution to the problem, and so we have 
followed the indexing used in [8]. 

2. Reference jSJ contains several hundred sequences related to the tooth- 
pick problem, many more than could be mentioned here. For a full list see 
[13] and also the entries in the index to [5] under "cellular automata." 

3. The computer language Matheinatica@ has a collection of com- 
mands that can often be used to display structures produced by cellular 
automata and to count the ON states. For example, the command 

Map [Function [Apply [Plus .Flatten [#1] ] ] , Cellular Automaton [{ 
686,{2,{{0,2,0},{2,1,2},{0,2,0}}},{1,1}},{{{1}},0},200]] 
produces the first 200 terms of the sequence U{n) giving the number of ON 
states in the Ulam-Warburton cellular automaton discussed in S|6] 
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Figure 4: Stages through 7 of the evolution of the corner structure. 
The numbers of toothpicks in the successive stages, C(0), . . . , C(7), are 
0,1,3,6,9,13,20,28. 



2. The corner sequence 

In order to understand the toothpick structure, it is helpful to first consider 
what happens if one quadrant of the plane is excluded. We impose the rule 
that no toothpick may cross into the third quadrant of the plane, and only 
ends of toothpicks may touch the negative x- or y-axes. At stage 0, we 
place a half-toothpick extending horizontally from the origin to the point 
(|, 0). The structure is then allowed to grow using the rule for the original 
toothpick sequence. The corner sequence is relevant because it describes 
how the main toothpick structure grows. 
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Let c{n) {n > 1) denote the number of toothpicks added at the nth stage, 
with c(0) = 0, and let C{n) := J27=o total number of toothpicks 

after n stages. These are respectively entries A152980 and A153006 in [S]. 

Fig. |4] shows stages through 7 of the evolution of the corner structure. 
Note that the first toothpick added, at stage 1 (with midpoint at the end 
of the initial half-toothpick), matches the initial toothpick of the original 
toothpick sequence, except that it is shifted a half-unit to the left. The 
initial values of c(n) and C{n) are shown in Table pj 
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Table 2: The corner sequences c{n) and C(n) for < n < 39. 






Figure 5: The corner toothpick structure after 2^ — 1 stages, for /c = 2. 3, 4. 



An examination of Fig. [4] and pictures of later stages in the evolution 
reveals that after 2^ — 1 stages (for A; > 2) the structure consists of an 
essentially solid rectangle of toothpicks with one quadrant removed. The 
first few cases are shown in Figs. [5j More precisely, we have: 
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Theorem 1. After 2^ — 1 stages, for k > 2, the corner toothpick structure 
is bounded by a rectangle of dimensions {height x width) — {2'^~^ — |) x 
^2^-1 _ X) with the lower left (2'^~^ ~ |) ^ (2*^"^ — 1) corner removed and 
with an additional half-toothpick protruding downwards from the lower right 
corner, in which all the boundary edges are solid rows of toothpicks except 
for the top edge which contains no horizontal toothpicks, with a row of 2^^^ 
exposed vertical toothpick ends along the top edge, and with no exposed 
toothpick ends in the interior. Furthermore, for k > 2, the number of 
toothpicks added at the successive stages while going from stage 2^ to stage 
2^=+! — \ is given by: 

[2^-1 + 1, if 1 = 0; 

c{2'' +i) = l2c{i}+c{i + l), ifi = \,...,2^-2; (2) 

\2c{i) + c{i + l)-l, if i = 2^-1. 

Proof. We use induction on k. The case fc = 2 is readily checked (cf. 
Figs. |4j [5]). Suppose the theorem is true for k, so that after the first 2*"' — 1 
stages we have the structure described in the theorem. We consider the 
next 2'' stages in the evolution of the bottom right quadrant and the top 
two quadrants separately. 

First, the bottom right quadrant looks like the starting configuration 
for the corner structure, with its protruding half-toothpick, except rotated 
clockwise by 90°. So by the induction hypothesis, after further 2*^ — 1 stages 
we reach a 90°-rotated copy of the (2^^ — l)-stage structure. One further 
step then fills in the right-hand edge, leaving a half-toothpick protruding 
downwards from the bottom right corner. Second, consider what happens 
to the top half of the structure. At the first step, the 2*"' — 1 vertical exposed 
toothpick ends will be covered, producing overhanging half-toothpicks at 
the left- and right-hand ends of the top edge. Again these look like the 
starting configuration for the corner structure, with the first quadrant a 
mirror image of the second quadrant. So again, by induction, after a further 
2*^ — 1 steps we reach the top half of the desired structure for k + 1. This 
completes the proof of the first assertion of the theorem. (The process is 
depicted schematically in Figs. [6] [t] and [s] ) 

The recurrence formula ([2]) now follows by keeping track of the number 
of toothpicks that are added at successive steps as we progress from stage 
2'= to stage 2''+^ - 1. ■ 

It is worth remarking that this growth in three quadrants, one of which 
is a step ahead of the other two, is responsible for the terms of the form 

2fiz) + fiz + 1) (3) 

which appear in the recurrences in Theorems [Ij [2] and |4] 
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Figure 6: Schematic of corner structure after 2*^ — 1 stages. 
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Figure 8: Schematic of corner structure after 2*^+^ — 1 stages. 




Figure 9: Schematic of evolution of toothpick structure: after 2^ stages 
(left-hand figure), after a further 2*^ — 1 stages (center) and after one more 
stage (right-hand figure). 
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3. The toothpick sequences 



Similar recurrences hold for the toothpick sequences t{n) and T[n). 

Theorem 2. For the toothpick structure discussed in f|7| the number of 
toothpicks added at the nth stage is given by t(0) = 0, t(l) — 1, and, for 
k > I, 

2t{i)+t{i + l), i/i = l,...,2'=-l. 



ti^'+^) = {L:. . . I '. ok . (4) 



Proof. An inductive argument similar to that used in the proof of Theorem 
[ijshows that after 2*"' steps, for k > 2, the toothpick structure is bounded by 
a 2*^^^ X (2*^^^ — 1) rectangle, with half-toothpicks protruding horizontally 
from the four corners, in which all the boundary edges are solid rows of 
toothpicks, and with no exposed toothpick ends in the interior. (The cases 
k — 2 and 3 can be seen in Fig. [l]) In the induction step, each quadrant 
grows like a suitably rotated version of the corner structure. (The evolution 
is depicted schematically in Fig. |9]) The recurrence formula Q now follows 
by keeping track of the number of toothpicks that are added as we progress 
from stage 2'' to stage 2'"'+^ — 1. ■ 

Corollary 3. For the toothpick sequence T{n), we have T'(O) = 0, and, for 
k>0, 



\T{2'')+2T{i)+T{i + l)-l, ifi = l, ' 
Proof. This follows easily from T{n) = X]"=o^(*)' 



J2 t(i) = 2'=(2'=+^ - 1) 

1=2''- 



and Q. We omit the details. ■ 

A convenient way to visualize the recurrences (|2|, (|4| and (§ is to write 
the sequences c(n), t{n) and T(n) as triangular arrays, with 1, 1, 2, 4, 8, 
16, 32, . . . terms in the successive rows. For example, the initial terms of 
the t{n) sequence are shown in the array in Table [sl The row labeled 8, for 
instance, begins with t{8) — 8, and then, using H) and referring back to 
the top of the triangle, continues with the values t{9) = 2t{l) + t{2) = 4, 
t(10) = 2i(2) + i(3) = 8, t(ll) = 2t(3) + <(4) = 12, and so on (a kind of 
"bootstrap" process). 

To see a direct connection between the toothpick sequences and the 
corner sequences, it is convenient to define Q{n) := (T(n) — 3)/4 (n > 3), 
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12 


16 28 


32 


16 


16 4 8 


12 


12 


16 28 


32 20 16 28 36 40 60 88 



Table 3: Initial terms of toothpick sequence t[n) arranged in triangular 
form. 

with Q(0) — (5(1) = 0. This is the number of toothpicks whose centers 
are in the interior of the first (or second, third or fourth) quadrants of the 
toothpick structure. Also let q{n) := Q{n) — Q{n—\) (n > 1) with q{0) ~ 0. 
The argument used in the proof of Theorem [T] shows that 

C{n) = 2Q{n) +Q{n+l) + 2, for n > 2 . (6) 

Hence by taking differences we have q{n) — t(n)/4, 

c{n) ^ 2q{n) + q{n + l), for n > 3 , (7) 

and 

c(n) = -i(n) + -<(n + 1), for ?i > 1 . (8) 
4. Rectangles in the toothpick structure 

Examination of Figs. [T||3] suggests that, after any finite number of stages, 
the toothpick structure divides the plane into an unbounded region and a 
number of squares and rectangles (and no other closed polygonal regions 
appear). Let R{n) denote the number of squares and rectangles in the 
toothpick structure after n stages, and let r(n) := R{n) — R{n — 1) be 
the number of squares and rectangles that are added at the nth stage. 
Similarly, let p{n) be the number of squares and rectangles that are added 
to the corner structure at the nth stage. The initial values of p(n), r(n) 
and R{n) are shown in Tabic |4] (these are entries A168131, A160125 and 
A160124 in 0). 

Then an inductive argument, similar to that used to establish Theorem 
[T] shows the following. 
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n 
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10 


11 


12 


13 


14 


15 


p{n) 
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18 


19 


r{n) 











2 


2 
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10 


6 
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8 


4 


4 


20 


30 


R{n) 











2 


4 


4 


8 


18 


24 


24 


28 


36 


40 


44 


64 


94 



Table 4: The sequences p{n), r{n) and R{n) for < n < 15. 



Theorem 4. All internal regions in the corner and toothpick structures are 
squares and rectangles. Furthermore, p{0) = p{l) = 0, p{2) = 1, p{3) ~ 2, 
and, for k > 2, 

2'=-i-l, z/z = 0; 

o(2'+^)-^^P^'^+P^' + ^^' z/z = l,...,2'=-3; 

^'^-'^2p{^)+pi^ + l) + l, ^J^ = 2>^-2. 
.2p(i)+p(z + l) + 2, i/i = 2'^-l; 



an. 



d r(G) = r(l) = r(2) = 0, r(3) = 2, anrf, for k>2, 

(2^-2, ifi^G- 
r{2^ + i)=Up{i), z/z = l,...,2'=-2; (10) 

yAp{i) + 2, i/i = 2'=-i. 

We omit the proof. 

5. The fractal-like structure 

The recursive structure established in the proofs of Theorems [T] and [2] also 
explains the fractal-like appearance of the toothpick array. After applying 
one round of the corner recursion to each quadrant and then rescaling, we 



have the transformation shown schematically in Fig. 10 (an "F" is used 



to indicate orientation of the various pieces), and in a specific example in 



Fig. 11 Note that four of the blocks (the sideways "F"s) are shifted by 
a half-toothpick towards the center. Because of this shift the toothpick 
structure is not strictly self-similar (cf. [3]) and so is not a true fractal. 
The same is true for all the structures we will meet in this paper: they 
have a fractal-like growth, but are not strictly self-similar. 

A plot of T(0), . . . , T{2^ — 1) for increasing values of k shows that the 
points lie roughly on a parabola, with irregularities caused by the fractal- 



like behavior (see for example Fig. 12). Benoit Jubin [S] has investi- 
gated limT(n)/n^ and lim T(n)/n^. His numerical results suggest that 
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Figure 10: Fractal-like transformation recursion, general step. 



Figure 11: Fractal-like recursion going from stage 8 to stage 16. 
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Figure 12: Plot of toothpick sequence r(0), . . . , T(1023). 

limr(n)/n^ — 2/3, with local maxima at the values n — 2*^ — 1, and 
]hnT{n)/n'^ w 0.4513058 with local minima at the following values of n: 

1, 2, 5, 12, 21, 44, 89, 180, 362, 728, 1459, 2921, . . . 

(see A170927 for further terms). His upper limit can be established from 
Corollary |3j 

Theorem 5. For n > 1, 

^<^^, (11) 

with equality if and only ifn = 2^ — \ for some k. Hence limT(n)/n^ — 2/3. 

Proof. The result is true for n < 3, and for n > 3 for the special values 
n^2'' - 1, when T{n) ^ {2^ - l){2^+^ - l)/3, so 

T{n) ^2^1^ 
rfi 3 3?! ' 

and n = 2^, when T{n) = [2"^^+^ + l)/3, so 



T(n) 2 12 1 

——L = - -I < - H 

3 "irfi 3 3rt 
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For the general case we use induction, and assume that (11) holds for all 
n < 2*^. Let n = 2^= + i, 1 < i < 2*^ - 2. Then (fill) follows from ((sl) and the 



induction hypothesis. ■ 

Jubin also observes that there is a continuous function on [0, 1] that 
describes the asymptotic behavior of T{n). This is the function whose 
graph is the Hausdorff limit of the finite sets Ek consisting of the points 
(x = jkJ[x) = ^) for n = 2*= + < i < 2'=. This function takes the 
value I at a; = and x = 1, and has its minimum at around (0.427451, 
0.4513058). It is non-differentiable at the dyadic rational points between 
and 1. Figure [T3| shows i?i4. 



Figure 14: Stages 1 through 8 of the evolution of the Ulam-Warburton 
structure. The numbers of ON cells in the successive stages, C/(l), . . . , J7(8), 
are 1,5,9,21,25,37,49,85. 



6. The Ulam-Warburton cellular automaton 



As we will see in pO| the toothpick structure can be modeled by a cellular 
automaton on a planar graph. In this section, we consider a simpler example 
of the same type, the arrangement of square cells generated by the Ulam- 
Warburton cellular automaton (Ulam [17], Singmaster TS", Stanley and 
Chapman |15j, Wolfram ^2] p. 928]). The cells are the squares in an infinite 
square grid, and the neighbors of each cell are defined to be the four squares 
which share an edge with it. (This is the von Neumann neighborhood of 
the cell, in the notation of [7].) At stage 0, no cells are ON. At stage 1, 
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a single cell is turned ON. Thereafter, a cell is changed from OFF to ON at 
stage n if and only if exactly one of its four neighbors was ON at stage n — 1. 
Once a cell is ON it stays ON. This is "Rule 686" in the notation of [5], [22] • 
Let u{n) [n > 0) denote the number of cells that are changed from OFF 
to ON at the nth stage, and let U{n) :— J^^^o ^'^ ^^'^ total number of 
ON cells after n stages. The initial values of u(n) and U{n) are shown in 
Table [5) These sequences are respectively entries A147582 and A147562 in 
[5j. Fig. [l4| shows stages 1 through 8 of the evolution of the this structure. 



(As is suggested by Fig. 14 and more particularly by the movie linked to 



entry A147562, this structure also has a fractal-like growth.) 



n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


u{n) 





1 


4 


4 


12 


4 


12 


12 


36 


4 


U{n) 





1 


5 


9 


21 


25 


37 


49 


85 


89 


n 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


u{n) 


12 


12 


36 


12 


36 


36 


108 


4 


12 


12 


U{n) 


101 


113 


149 


161 


197 


233 


341 


345 


357 


369 


n 


20 


21 


22 


23 


24 


25 


26 


27 


28 


29 


u{n) 


36 


12 


36 


36 


108 


12 


36 


36 


108 


36 


U{n) 


405 


417 


453 


489 


597 


609 


645 


681 


789 


825 


n 


30 


31 


32 


33 


34 


35 


36 


37 


38 


39 


u{n) 


108 


108 


324 


4 


12 


12 


36 


12 


36 


36 


U{n) 


933 


1041 


1365 


1369 


1381 


1393 


1429 


1441 


1477 


1513 


n 


40 


41 


42 


43 


44 


45 


46 


47 


48 


49 


u(n) 


108 


12 


36 


36 


108 


36 


108 


108 


324 


12 


U{n) 


1621 


1633 


1669 


1705 


1813 


1849 


1957 


2065 


2389 


2401 



Table 5: The sequences u{n) and U{n) from the Ulam-Warburton cellular 
automaton, for < ?i < 49. 



Theorem 6. (i) The number of cells that turn from OFF to ON at stage n of 
the Ulam-Warburton cellular automaton satisfies the recurrence u(0) = 1, 
u(l) = 1, and, for k > 0, 

[3u{i), tfi = l,...,2'^-l. 
(ii) There is an explicit formula: u(0) — 0, u(l) ~ 1 and 

= 4-3"*("~i) - 1, n> 2, (13) 
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where wt(n), the "binary weight" of n, is the number ofl's in the binary 
expansion of n (entry A000120 in 
(Hi) The u[n) have generating function 

^(-^ + ^11(1 + 3^''))- (14) 

fe>0 

Proof. Part (i) follows by an inductive argument similar to that used 
in the proofs of Theorems [T] and [2j The appropriate "corner sequence" is 
A048883, in which the nth term is 3wt("-i) (n > l), with partial sums given 
by A130665. Part (ii) follows from (i) by induction on n. Part (iii) follows 
from the generating function for A048883, which is nfe>o(-'^ ^a;^*^). ■ 



Remarks. 1. Now the corner sequence has three quadrants that grow in 
synchronism, so the 2c(z) + c{i + 1) terms in ^ are replaced by the 3u(i) 
term in ( 12 1. 

2. Parts (i) and (ii) of the theorem can be found in Singmaster [13] and 
Stanley and Chapman [13], and part (i) at least was probably known to J. 
C. Holladay and Ulam. On page 216 of [T7], Ulam remarks that for certain 
structures similar to this one (exactly which ones is left unspecified), J. C. 
Holladay showed that "at generations whose index number n is of the form 
n = 2*^, the growth is cut off everywhere except on the 'stems', i.e. the 
straight lines issuing from the original point." This is certainly consistent 
with the recurrence (12). 

3. Two properties of the Ulam-Warburton structure given in 15] are 
worth mentioning here, (i) When considered as a subgraph of the infinite 
square grid, the structure is a tree. This is also true for the toothpick 
structure — see \ 10 (ii) Let n — 1 = Yll'=i 2''" (ri > r2 > • • • > r^, > 0) be 



the binary expansion of n — 1. Then a necessary and sufficient condition 
for the cell at P = (x, y) G Z x Z to be turned from OFF to ON at stage 
n > 1 is that P = Er=i2'''u,, where G {(-1, 0), (1, 0), (0, -1), (0, 1)}, 
subject to Vi 7^ — Ui-i for i > 1. We have no such characterization of the 
toothpicks added at the nth stage. It is a consequence of this (although it 
is not mentioned in [121) ^^lat the cells that are turned ON at some stage 
are the cells (cc, y) with a; = or ?/ = 0, and the cells with xy for which 
the highest power of 2 dividing x is different from the highest power of 2 
dividing y. Again we know of no analog for the toothpick structure. 



7. Leftist toothpicks 

Stimulated by Theorem [6j we set out to look for analogues of (131 and 
for the toothpick sequence t{u). Our first attempt was a failure, but led to 
an interesting connection with Sierpinski's triangle. 
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We define the "leftist toothpick" structure as follows. We start with 
a single horizontal toothpick at stage 1, and extend the structure using 
the toothpick rule of |T] except that if a toothpick is horizontal, a new 
toothpick can be added only at its left-hand end. Let l{n) {n > 1) denote 
the number of toothpicks added at the nth stage, with 1{Q) = 0, and let 
L{n) := X]"=o^(*) total number of toothpicks after n stages. These 

are respectively entries A151565 and A151566 in ^8,. The initial values of 
l{n) and L{n) are shown in Table [6j Figure 15 shows the first 15 stages 
of the evolution (the starting toothpick is the apex of the triangle, at the 
right). 



n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


l{n) 





1 


1 


2 


2 


2 


2 


4 


4 


2 


2 


4 


4 


4 


4 


8 


L{n) 





1 


2 


4 


6 


8 


10 


14 


18 


20 


22 


26 


30 


34 


38 


46 



Table 6: The leftist toothpick sequences l{n) and L{n) for < n < 15. 



Figure 15: Stages 1 through 15 of the evolution of the leftist toothpick 
structure. 

The reason for investigating this structure is that, at least for the early 
stages, the part of the toothpick structure of f|T] in the 90° sector x — 
2 < y < 2 — X is essentially equal to the leftist structure. This breaks 
down, however, at stage 14. Nevertheless, the leftist structure has some 
interest. For if we rotate the structure anticlockwise by 90° and erase all the 
horizontal toothpicks, we obtain a triangle in which the vertical toothpicks 
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correspond exactly to the positions of the Is in Sierpihski's triangle (i.e., 
Pascal's triangle read modulo 2 [4], [10], [12], [22l Chap. 3]), and the gaps 
between the vertical toothpicks to the Os. Once observed, this is easy to 
prove. 

Gould's sequence, entry A001316 in [5], gives the number of odd entries 
in row n of Pascal's triangle, which is 2'^*'^"^ with generating function 
]^^>p(l + 2x^ ). Allowing for the different offset, we conclude that the 
leftist toothpick sequence l{n) is given by l{2n - 1) = /(2n) = 2**("-i). 



8. Generating functions 

The following theorem suggests why generating functions of the form ([T]) 
arise in connection with recurrences of the form ([2]), Q. 

Theorem 7. Given integers a, /3, 7, 5, let the Taylor series expansion of 

x{a + /3x)Y[{l + jx^'°-'^ +Sx^'') (15) 
fc>i 

be a(0) + a(l).T + a{2)x'^ + • • • . Then we have a(0) = 0, a(l) = a, and for 
n>2, 

(aj + lSd^-'^ ifi = 0; 

a{2'' +i) = <5a{i) +-fa{i + l), if i ^ 1, . . . ,2'' - 2; (16) 

[Sa{i) + ■yaii + 1) - aj^ , ifi = 2^-l. 

Proof. The generating function is 

x{a + (ix){l + 7a; + 5x^){l + ^x^ + 5x'^){l + ^x'' + bx^) ■■■ . 

Consider the coefficient of x^'^, say. The only way to build up x^'^ is to 
combine the term 6x^ with a(2)a;^, or the term 7a;'' with a{3)x^. Hence 
a(2'^ + 2) = da{2) + 70(8). Similar arguments holds for the general case, 
although adjustments are needed when i = or 2*^ — 1. We omit the 
details. ■ 

An analogous result holds if the product in (151 starts at A: = 0. 
Notice the resemblance between Equations ([2|), Q and (16 1. In partic- 
ular, we have: 

Theorem 8. The generating function for the corner sequence c{n) is 

C{x) := a;(l + a;) J|(l + x^"-'^ + 2x^'') . (17) 
fc>i 
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The generating function for the toothpick sequence t{n) is 

r{x) := ^{1 + Ax{l + X) n (1 + + -ix^")} 

^ k>l 

= ^{l + 2xl[{l + x''-'+2x'')}, (18) 

^ i;>0 

and therefore the generating function for the toothpick sequence T{n) is 

(i + 2.na + ^''-' + 2.^-)). (1.,) 

Proof. For the first assertion, we set a = /3 = 7 = 1, ^ = 2 in Theorem 
[7] and use ([2]). For the second assertion, we note that ([t]) imphes that the 
generating functions for c(n) and q{n) are related by 

C{x)^x + x^ + (2 + -)Q{x), (20) 
a; 

and by definition we have 

T{x)=x + 2x'^ + Q{x). (21) 



Ehminating Q{x), we obtain (18). 



Remark. Equation (19 1 was conjectured by Gary W. Adamson [T]. Con- 
sider the sequence 1, 1, 2, 1, 3, 4, 4, 1, 3, 4, 5, . . . with generating function 

l[{l + x^'-'+2x''') 

k>l 

(entry A151550 in [S'). Adamson discovered that if this sequence is con- 
volved with the sequence 1, 2, 2, 2, 2, 2, . . ., the result appeared to coincide 
with the corner sequence C(n). When expressed in terms of generating 
functions, his conjecture is essentially equivalent to (19). 

9. Explicit formulas 

A second comment in |5], this time from Hagen von Eitzen, was instru- 
mental in the discovery of explicit formulas for many of these sequences. 
Von Eitzen [18] contributed the sequence 2, 3, 3, 3, 5, 6, 4, 3, 5, 6, 6, . . . with 
generating function nfe>o(-'^ ~^ ) to j8] (it is entry A160573) and 
provided an elegant explicit formula for the nth term: 



m>0 



E r*^i^"'^)- (22) 
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This can be generalized. For this it is convenient to omit the initial 
linear factors from ([ij but to start the product at A: = 0. 

Theorem 9. Let the Taylor series expansion of 

\{{l+^x^"-^ + 5x^'') (23) 

A->0 

be a(0) + a[l)x + a{2)x'^ + • • • . Then 

a{n) = ^-jwt(„+™)-™ + ^)\ (24) 

m>0 \ m, J 



Proof. (Based on von Eitzen's proof of (22).) First, observe that 



[](l + dV')-^J-*(")a;", (25) 



fc>0 n=0 



since when getting a term x", we pick up a factor of 5 for every 1 in the 
binary expansion of n. When we expand 

Wil+^x^'-'+Sx^"), (26) 

fc>0 



instead of the product in (25 1, each time we replace a term 5x by 7X , 
we lose a factor of x in the product, but we gain because there may be 
several ways to choose the factors in which to do the replacement. Suppose 
we do this replacement in m of the terms in ( [26| . Then we must increase 
n to n + m, we gain by a factor of (**^"^™''), but we have to replace m 
factors of 5 by 7s, for a net contribution of ^'"5wt(n+m)-m^wt(n+m)^ ^-^^ 
sum. ■ 

Remark. Note that there are only finitely many nonzero terms in the 



summations (22 1 and (24). For large n the number of nonzero terms is 
roughly log2 n. More precisely, the number of nonzero terms for any n is 
given by entry AIOO66I in [5]. 

We can use Theorem [9] to obtain an explicit formula for the toothpick 
sequence t{n). 

Theorem 10. 

V ' Ufc+i^ ifi = 2^-l. ^ ' 
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Proof. This theorem is an instance where our decision to start enumera- 
tions at n = rather than n = 1 (as discussed in Remark 1 in SjT]) causes 
comphcations. The result would be more elegant if we had made the other 
choice. So, just for this proof, let us define i{n) — t{n + 1) for n > 0. Then 
from Theorem [2j i{n) satisfies the recurrence 

f(2^+.) = ^;2 + ^"(^ + l)' if-0;----2'^-2; (28) 
V ' l2'^+\ ifi = 2'=-l. ^ ' 



k 


terms 2'=,2'= + 


1,.. 




n _ 


1 
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4 8 


12 8 












8 


4 8 


12 12 


16 


28 


32 


16 




16 


4 8 


12 12 


16 


28 


32 


20 16 28 36 40 60 8 


8 80 



Table 7: Initial terms of sequence t{n) arranged in triangular form. 

The initial terms of the t{n) sequence are shown in triangular array form 
in Table [tI The initial 2*^ — 1 terms of the fcth row are the initial 2*^ — 1 
terms of the next row, so the rows converge to a sequence 

4, 8, 12, 12, 16, 28, 32, 20, 16, 28, 36, . . . 

that we will denote by F{n), n > (this is entry A147646 in [8]). It 
follows from the definition that F{n) is defined by the recurrence F{0) = 4, 
F(l) = 8, F(2) = F(3) = 12, and for n > A, 

(2F{i) + F{i + l), if0<i< 2*^-3; 
F{2'' + i)= l2F{i)+F{i + l)~A, Hi = 2'' -2; (29) 
[2'=+2+4, if i = 2*^-1. 

Also, for < i < 2*^ - 2, 

<(2'= +i)= F{i) . (30) 

Comparing (29 1 with Theorem [7j we see (taking a = /3 = 4, 7 = 1, 5 = 2 
in that theorem) that F{n) has generating function 

4(1 + x)ll{l + x^"-' + 2x^") = 2 + 2^''"' + 22:^') • (31) 

k>l k>0 
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It follows from Theorem |9] that 

= 2 J2 2-'(^+")-" ("^^^'^ '"^ ) . (32) 



m>0 



Then (30) and (32 1 imply (271 



10. Cellular automata defined on graphs 

Cellular automata defined on graphs were introduced by von Neumann and 
Ulam in 1949 [16], and the toothpick structure {^^, the corner toothpick 
structure the leftist structure (Q, the Ulam-Warburton cellular au- 
tomaton (0S, etc., can all be described in this language. Let G be an 
infinite directed (or undirected) graph with finite in-degree and out-degree 
(or degree) at each node. Nodes are either in the OFF state or the ON state, 
and we just need to give a rule for deciding when the nodes change state. 

For the toothpick structure, the nodes of G are the vertices {x,y) E Z x Z 
of the square grid. There are two kinds of nodes: "even" nodes, with x + y 
even, which have edges directed to nodes {x,y± 1), and "odd" nodes, with 
x + y odd, which have edges directed to nodes {x ± 1, y). Initially all nodes 
are OFF, at stage 1 we turn ON node (0,0), and thereafter a node turns ON 
if it has an incoming edge from exactly one ON node. This is easily seen to 
be equivalent to the toothpick structure, with the even (resp. odd) nodes 
representing the midpoints of vertical (resp. horizontal) toothpicks. The 
induced subgraph of G joining the ON nodes is a directed tree (and remains 
a tree if the arrows on the edges are removed) . 

For the Ulam-Warburton cellular automaton, of course, G is the undi- 
rected graph with vertices {x, y) S Z x Z, with each node connected to its 



four neighbors (or, in the case of the cellular automaton analyzed in { 14 
its eight neighbors). As already mentioned in f|6j the induced subgraph 
joining the ON nodes is also a tree. 

For the natural generalization of the Ulam-Warburton cellular automa- 
ton to higher dimensions, with G = Z'^, 5 > 1, and each node adjacent to 
its 2d neighbors, there is an analog of Theorem |6] The number of cells that 
turn from OFF to ON at stage n is given by u(0) — 0, u{l) = 1, and, for 
n > 2, 

w(n) = 2d(2d- l)**("^i)-i (33) 

(Stanley and Chapman [15]; see also entries A151779, A151781 in |8]). On 
the other hand, for the face-centered cubic lattice graph, where each node 
has 12 neighbors, there is no obvious formula (see A151776, A151777). 

The number of possibilities is endless: see [Hj and [55] for further ex- 
amples. 
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Figure 16: Stages 1 through 4 of the evolution of the T-toothpick structure. 
11. T-shaped toothpicks 

One may consider "toothpicks" of many other shapes. Here we consider just 
one example: we define a "T-toothpick" to consist of three line segments 
of length 1, forming a "T". Using an obvious terminology, we will speak of 
the "crossbar" (which has length 2) and the "stem" (of length 1) of the T, 
and refer to the midpoint of the crossbar as the "midpoint." A T-toothpick 
has three endpoints, and an endpoint is said to be "exposed" if it is not the 
midpoint or endpoint of any other T-toothpick. 

We start at stage 1 with a single T-toothpick in the plane, with its 
stem vertical and pointing downwards. Thereafter, at stage n we place a 
T-toothpick at every exposed endpoint, with the midpoint of the new T- 
toothpick touching the exposed endpoint and with its stem pointing away 
from the existing T-toothpick. Figure [16] shows the first four stages of 
the evolution. Let T{n) denote the number of T-toothpicks added to the 
structure at the nth stage (this is A160173). 

Theorem 11. We have r(0) = 0, t(1) = 1, t(2) ^ 3, and, for n>3, 

r(n) = ^{3wt("-i) _^ 3wt(n-2)| _^ ^ ^34^ 
3 

Proof. This is easily established from Theorem [6] by observing that the 
structures in the four quadrants defined by the the initial T are essentially 
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equivalent to copies of the Ulam-Warburton structure. The copies in the 
first and second quadrants are one stage behind those in the third and 
fourth quadrants. ■ 

The analogous structure using Y-shaped toothpicks has resisted our 
attempts to analyze it — see A160120 in [5]. 

12. The "Maltese cross" or Holladay-Ulam structure 

On page 217 of fTT!, Ulam discusses another structure that he and J. C. 
HoUaday had studied. To construct this, one first builds the infinite Ulam- 
Warburton structure described in f|6j and then replaces each ON cell by a 
Maltese cross consisting of a central ON cell surrounded by four other ON 
cells. Now label the cells of the new infinite structure, starting by labeling 
the central square 1, then the four adjacent cells 2, and so on, always moving 
outwards from the center. Figure [T7| shows the cells with labels 1 through 
5. 
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Figure 17: Stages 1 through 5 of the evolution of the Maltese cross structure. 

HoUaday and Ulam then give a set of rules for a cellular automaton that 
will build up this structure, starting with the cell labeled 1. We found their 
rules (stated on pages 216 and 222) somewhat hard to understand, so it 
may be helpful to the reader if we give our version of them here. 

For a candidate cell to be turned ON it is necessary (though not sufficient) 
that it share exactly one edge with an ON cell of the previous generation. 
The two vertices of the candidate cell that touch that edge we will call its 
inner vertices, and the other two vertices we call its outer vertices. 
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The rules are as follows. Cells are in one of three states, OFF, ON or 
DEAD. Once a cell is ON or DEAD it remains in that state. Initially all cells 
are OFF, and at stage 1 a single cell is turned ON. 

At stage n, consider OFF cells X (the "candidates") that share at least 
one edge with an ON cell of the previous generation. If a candidate X is 
edge-adjacent to two ON cells, X is declared DEAD. If two candidates X and 
X' share an outer vertex, both X and X' are declared DEAD. If a candidate 
X is edge-adjacent to a DEAD cell, then X is declared DEAD, except that 
when n = 2 (mod 3), this rule obtains only if X is edge-adjacent to a cell 
that was declared DEAD at the previous generation. If a candidate is not 
excluded by these conditions, it is turned ON. 

Figure [18] is an annotated picture showing the first eight stages of the 
evolution of this structure in the first quadrant, with letters identifying the 
first few DEAD cells. The letters 'a' and 'd' indicate candidate cells that are 
DEAD because they are edge-adjacent to two ON cells, and 'b' and 'e' indicate 
candidate cells that are DEAD because they are edge-adjacent to a DEAD cell 
and n is not congruent to 2 (mod 3). The '5' cells not on the main axes are 
ON because 5 = 2 (mod 3). The two 'c' cells are DEAD because they would 
share a common outer vertex; and so on. 
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Figure 18: Annotated version of stages 1 through 8 of the Maltese cross 
structure in the first quadrant. 

Of course, the original definition of this structure makes it easy to give 
a formula for the number of cells, m{n), say, that are added at the nth 
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stage. From ([Is]) we have m{0) — 0, m{l) — 1, m{2) = 4, and, for t > 1, 

m{3t) = m(3i + 1) = 4-3**(*)-\ m{M + 2) = 4-3^'(*). (35) 

(This is entry A151906 in "W.) 

We conclude this section by hsting some related cellular automata stud- 
ied by Ulam and his colleagues HoUaday and Schrandt in [TT] and [T7] . Since 
we have not even found recurrences for them, we will give no details. 

Reference [llj mentions a cellular automaton which is intermediate in 
complexity between the Ulam-Warburton structure and the Maltese cross 
structure discussed above. This Schrandt-Ulam cellular automaton is de- 
scribed in entries A170896 and A170897 in [8j; another version is given in 
A151895, A151896. 

Analogous structures based on triangles or hexagons rather than squares 
are described in Examples 3 through 6 of [17] (see also Wolfram ^22, p. 371]), 
and in entries A151723, A151724, A161644, A161645 of [5]. We invite the 
reader to find recurrences or generating functions for any of them. 



13. Rule 942 



On page 928 of 22 , Wolfram considers (among other examples) what hap- 
pens if the rule for the Ulam-Warburton cellular automaton is modified so 
that a cell turns ON if and only if either exactly one or all four of its four 
neighbors is ON (this is "Rule 942" in the notation of [9], [22]). 

Let w{n) denote the number of cells that are changed from OFF to ON 
at stage n. Since the four-neighbor part of the rule is invoked only after 
an OFF cell is completely surrounded by ON cells, w(n) > u{n) for all n. 
In fact, w' {n) :~ w{n) — u(n) is always a multiple of 4, and w^n) — u{n) 
except when n = 1 (mod n). Table [s] shows the initial values of w{n)^ u{n), 
w'{n) and 6{n) :== |(w(4r7. + 1) - u{4n + 1)) (cf. A169648 and A169689 in 
[Si). 
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11 


12 


12 


14 


15 


w(n) 
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4 


12 
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12 


12 


36 


28 


12 


12 


36 


28 


36 


36 


u(n) 
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12 


12 


36 
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12 


12 


36 


12 


36 


36 


w'{n) 

















4 











24 











16 








S{n) 





1 


6 


4 


24 


4 


20 


12 


84 


4 


20 


12 


76 


12 


60 


36 



Table 8: The sequences w{n), u{n), ■w'{nQ := w(n) — u{n), 6{n) for < n < 
15. 



There is a simple explicit formula for 6{n) and hence, via ( 13 1, for w{n 
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Theorem 12. We have 5{Q) = 0, 5{1) = 1, ^(2) = 6. For n > 3, let 

n^2^ + j with l<j< 1^, where j = 2'"(2; + 1) (say). Then 

S{n) = 4(3™+i - 2"+i)3*'(') , (36) 
except that if j = 2^ then 

8{n) = 4-3'=+i - 3-2'^'+i . (37) 
We omit the proof, which is similar to those of Theorems [l] and [2] 

14. Square grid with eight neighbors 

Our final example is also based on the Ulam-Warburton cellular automaton, 
except that now we take the neighbors of a cell to consist of the eight cells 
surrounding it. (This is the Moore neighborhood of the cell, in the notation 
of (Z)-) Otherwise the rules are the same as in ^joj a cell turns ON if exactly 
one of its eight neighbors is ON. 

Let v{n) (n > 0) denote the number of cells that are changed from 
OFF to ON at the nth stage of the evolution, and let V{n) := X)"=o ^(*) ^® 
the total number of ON cells after n stages. The initial values of v{n) and 
V{n) are shown in Table ^ below. These sequences are respectively entries 
A151726 and A151725 iii[8]. Figure [l9] shows stages 1 through 8 of the 



evolution of the this structure. 

Since each cell now has two kinds of neighbors, it is perhaps not surpris- 
ing that this problem is more difficult to analyze than the Ulam-Warburton 
structure. In order to understand the growth, it is convenient to define two 
versions of "corner sequences," analogous to that introduced in f|2j 

So that we can refer to individual cells, we will label each square cell 
by the grid point at its upper left corner. That is, we define cell {i,j) to 
consist of the square {{x,y) ^ M. x M. \ i < x < i + 1, j — I < y < j}- 

For the first corner sequence, we exclude the third quadrant of the 
plane, and at stage 1 we turn ON the cell immediately to the right of that 



quadrant (see Fig. 20). More precisely, at stage 1, we turn ON the cell (0,0), 
and thereafter extend the structure using the eight-neighbor rule, with the 
proviso that after the first stage, no ON cell may be adjacent to any of the 
third-quarter cells — meaning the cells («, j) S Z x Z with i < — l,j < 0. 

The second corner sequence is similar to the first, except that at stage 1 
we turn ON the cell (0, 1), just up and to the right of the excluded quadrant 



(Fig. 21 1 



Let Vi{n) (resp. V2{n)) denote the number of cells that are changed 
from OFF to ON at the nth stage of the evolution of the first (resp. second) 
corner sequence. The initial values of ui(n) and f2(n) are also shown in 
Table [9] These sequences are respectively entries A151747 and A151728 in 
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Figure 19: Stages 1 through 8 of the evohition of the eight-neighbor struc- 
ture. The numbers of ON ceUs in the successive stages, V{1), . . . , V{8), are 
1,9,13,33,37, 57, 77,121. 



[8]. Figures 20 and 21 shows stages 1 through 5 of the evolution of the two 
corner sequences. 

The fohowing theorem gives recurrences for all three of these sequences. 

Theorem 13. The eight-neighbor sequences vi{n), V2(n) and v{n) satisfy 
the following recurrences: 

vi{0) = 0, wi(l) = 1, vi{2) = 3, i;i(3) = 5, and, for k>2, 



vi{2''+i) 



■(3fc + l)2'^-2_^l, ifi = 0; 

3-2'=-i+3, ifi = l; 

2vi(i)+wi(i + l), «/i-2,. 



(38) 



2vi{i) + vi{i + l)-l, ifi = 2'' 

V2{0) = 0, ^2(1) = 1, and, for k > 1, 

'3-2'^ -1, ifi = 0; 

V2{2''+i) = {v2{i) + 2vi{i + l), j/« = l,. 



,2'= - 2: 



(39) 



W2(») + 2ui(i + 1) - 2, ifi = 2^-l\ 
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44 


28 


60 


76 
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60 


84 


116 


116 
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581 
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897 
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41 


40 


61 
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29 


49 


61 


53 


29 


51 


71 


87 


107 



Table 9: The 8-neighbor sequences v{n) and V{n), and the two "corner 
sequences vi{n), V2{n), for < n < 23. 



v{0) = 0, v{l) = 1, and, for k>l, 

.ot- X [6-2'= -4, ifi = 0; 

[4w2(i), z/i = 1,...,2'= - 1. 

Again we omit the proof. We have not found generating functions or 
explicit formulas for any of these sequences. 



15. Acknowledgments 

We thank Benoit Jubin for telling us about his investigations into the 
asymptotic behavior of T{n) that were discussed in qSj and Gary Adam- 
son and Hagen von Eitzen for their contributions to Ts] which were men- 
tioned in S|8] and ^ We also thank Maximilian Easier, John Layman and 
Richard Mathar, who have made many contributions to [5] (especially new 
sequences or extensions of existing sequences) related to the subject of this 
paper. Finally, we thank Laurinda Alcorn for locating a copy of |11] . 



References 



[1] G. W. Adamson, Comment on entry A151550 in 'B], May 25, 2009. 



34 



[2] E. R. Berlekamp, J. H. Conway and R. K. Guy, Winning Ways for 
Your Mathematical Plays, A. K. Peters, Wellesley, MA, 2nd. ed., 4 
vols., 2001-2004. 

[3] K. Falconer, Fractal Geometry, Wiley, NY, 1990. 

[4] S. R. Finch, Mathematical Constants, Cambridge, 2003, p. 145. 

[5] H. Gingold, H. W. Gould and M. E. Mays, Power product expansions. 
Until. Math., 34 (1988), 143-161. 

[6] B. Jubin, Personal communication, January, 2010. 

[7] M. Mitchell, Computation in cellular automata: a selected review, pp. 
95-140 of T. Gramfi et al.. Nonstandard Computation, Wiley- VCH, 
Weinheim, 1996. 

[8] The OEIS Foundation Inc., The On-Line Encyclopedia of Integer Se- 
quences, |http : // oeis . org , 2010. 

[9] N. H. Packard and S. Wolfram, Two-dimensional cellular automata, J. 
Statistical Physics, 38 (1985), 901-946. 

[10] H.-O. Peitgen, H. Jiirgens and D. Saupe, Chaos and Fractals, Springer- 
Verlag, NY, 1992, pp. 408-409. 

[11] R. G. Schrandt and S. M. Ulam, On recursively defined ge- 
ometric objects and patterns of growth, Los Alamos Scien- 
tific Laboratory, Los Alamos, NM, Report LA-3762, Aug 16 
1967; published in A. W. Burks, editor. Essays on Cellular Au- 
tomata, Univ. 111. Press, 1970, pp. 238ff. Available on-line from 
|http r77library . lanl . gov/cgi-bin/getf ile?0035¥637 .pd?] 

[12] M. Schroeder, Fractals, Choas, Power Laws, W. H. Freeman, NY, 1991, 
Chap. 17. 

[13] D. Singmaster, On the cellular automaton of Ulam and Warburton, 
M500 Magazine of The Open University, #195 (December 2003), 2-7. 

[14] N. J. A. Sloane, Catalog of Toothpick and Cellular Au- 
tomata Sequences in the OEIS, published electronically at 
|http: //www. research, att . com/^njas/sequences/toothlist .html] 
2010. 

[15] R. P. Stanley (proposer) and R. J. Chapman (solver), A tree in the 
integer lattice. Problem 10360, Amer. Math. Monthly, 101 (1994), 76; 
105 (1998), 769-771. 



35 



[16] S. M. Ulam, Random processes and transformations, in Proc. Interna- 
tional Congress Mathematicians, Cambridge, MA, 1950, Amer. Math. 
Soc, Providence, RI, 1952, Vol. 2, pp. 264-275. 

[17] S. M. Ulam, On some mathematical problems connected with patterns 
of growth of figures, pp. 215-224 of R. E. Bellman, ed., Mathematical 
Problems in the Biological Sciences, Proc. Sympos. Applied Math., 
Vol. 14, Amer. Math. Soc, 1962. 

[18] H. von Eitzcn, Entry A160573 in [8], May 20, 2009. 

[19] E. W. Weisstein, Rule 30, published electronically at 
|http : //mathworld . wolfram . com/Rule30 . html (From MathWorld-A 
Wolfram Web Resource). 

[20] S. Wolfram, Statistical mechanics of cellular automata. Rev. Mod. 
Phys., 55 (1983), 601-644. 

[21] S. Wolfram, The Mathematica@ Book, Wolfram Media, Champaign, 
IL, 3rd. ed., 1996. 

[22] S. Wolfram, A New Kind of Science, Wolfram Media, Champaign, IL, 
2002. 



36 



